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Abstract 
We use simulated Hubble parameter data in the redshift range < z < 2 to explore the role 
and power of observational H{z) data in constraining cosmological parameters of the ACDM model. 
The error model of the simulated data is empirically constructed from available measurements and 
scales linearly as z increases. By comparing the median figures of merit calculated from simulated 
datasets with that of current type la supernova data, we find that as many as 64 further independent 
measurements of H{z) are needed to match the parameter constraining power of SNIa. If the error of 
H{z) could be lowered to 3%, the same number of future measurements would be needed, but then the 
redshift coverage would only be required to reach z — 1. We also show that accurate measurements 
of the Hubble constant Hq can be used as priors to increase the H{z) data's figure of merit. 

Subject headings: cosmological parameters — dark energy — distance scale — methods: statistical 
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1. INTRODUCTION 

The expansion of the Universe can be quantitatively 
studied using the results from a variety of cosmological 
observations, for example the mapping of t he cosmic mi- 
crowave background (CMB ) anisotropics (Spergel et al. 
2010), the measuremient ot baryon 



structed from the luminosity distances of SNIa ( Wang fc 



2007 



acous 



Komatsu et al. 



:ic oscillation (B aO) peaks ( Eisenstein et al.||2005 
Percival et al. 2 010), the linear growth ot large-scale 



structures (LSS) (|Wang fc Tegmark'2004), strong gravi- 



tational lensing ( \ang fc (Jhen_2009} , and measurements 
of "standard candles" such as the redshift-distance re- 




lationship _o£_t;5[£eJa^_supernovae (SNIa) (Riess et al. 

f^% ' _ ' ' ' 

hirianda 
observations is tne determmation ot tne Mubble param- 
eter H which is directly related to the expansion history 
of the Universe by its definition: H = a/ a, where a de- 
notes the cosmic scale factor and a is its rate of change 
with respect to the cosmic time. In practice, the Hubble 
parameter is usually measured as a function of the red- 
shift z, and the redshift is related to a by the formula 
a{t)/a{tf)) = 1/(1 -|- z) where t^ is the current time that 
is taken to be a constant. In the rest of this paper we 
will use the abbreviation "OHD" for observational H{z) 
data. 

Though not directly observable, H{z) can neverthe- 
less be deduced from various observational data, such 
as cosmological ages ("standard clocks") or sizes ("stan- 
dar d rods' ). The fo r mer m ethod has been illustrated 
by Jimenez fc Loeb (2002) and leads to OHD found 
from difierentiai ages ot galaxies. The latter meth od has 
been discu s sed by Blake & Glazcbrook 



2003) and|Seofc 



Eisenstein ( 2005[ ), which leads to H{z) found from BAO 
peaks. 

Moreover, the Hubble parameter as a quantitative 
measure of the cosmic expansion rate is closely related 
to cosmological distances. In particular, it can be recon- 
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Tegmark|[2005l IShafieloo et al.||2006)|Mignone fc Sarte 



mann|2008p "and GRBs ( [Liang et al.|2010 ). On the other 
hand, given the aforementioned availability of indepen- 
dent determination of H{z), and the expectation of more 
available data in the future, we are interested in the 
comparison of OHD and SNIa data in terms of their re- 
spective merits in constraining cosmological parameters. 
This interest of ours is expressed by the following two 
questions: Can future observational determinations of 
the Hubble parameter be used as a viable alternative to 
current SNIa data? If so, how many more datapoints are 
needed so that the cosmological parameter constraints 
obtained from H{z) data are as good as those obtained 
from SNIa distance-redshift relations? 

We attempt to illustrate possible answers to these two 
questions via an exploratory, statistical approach. This 
paper is organized as follows: we first briefiy summarize 
the current status of available OHD results in Section 
[2J Next, we show how the simulated H{z) datasets are 
used in our exploration in Section [3] and present the 
results from the simulated data in Section |4l In Section 
[5] we turn to the data expected from future observation 
programs. Finally, in Section [6] we discuss the limitation 
and implications of our results. 

2. AVAILABLE HUBBLE PARAMETER DATASETS 

Currently the amount of available H{z) data is scarce 
compa red with SNIa luminosity distance data 



Jimenez 



et al. (2003 ) first obtained one H{z) data point at z ~ 
0.1 from observatio ns of galaxy ages (henceforward the 
"JVTS03" dataset). [Simon et al.| ( [2005 J further obtained 
8 additional H{z) values up to z = 1.75 from the rel- 
ative ages of passively evolving galaxies (henceforward 
"SVJ05") and used it to const rain the redshift de pen- 



dency of dark energy potential. Stern et al. (2010a) ob 



tained an expanded dataset (henceforward SJVKSIO' 
and combined it with CMB data to constrain dark en- 
ergy parameters and the spatial curvature. Besides 
H{z) determinations from galaxy ages, observations of 
BAO peaks have also been used to extr act H{z) val- 
ues at low redshift ( Gaztanaga et al.[|2009 henceforward 
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Table 1 

Currently Available H{z) Datasets 



SVJ05 



SJVKSIO GCH09'' 



0.09'= 

0.17 

0.24; 

0.27 

0.4 

0.43 

0.48 

0.88 

0.9 

1.3 

1.43 

1.53 

1.75 



+0.06 
0.09 



+0.04 



69 ±12 
83 ± 8.3 

70 ±14 
87 ± 17.4 



117 ±23.4 

168 ± 13.4 
177 ± 14.2 
140 ± 14 
202 ±40.4 



69 ±12 

83 ±8 

77 ±14 
95 ±17 

97 ±62 
90 ±40 
117 ±23 
168 ±17 
177 ±18 
140 ± 14 
202 ± 40 



79, 



86, 



69 ±2.65 



45 ±3.1 



Note. — H{z) figures quoted in this tabic are 
in the units of km s~^ Mpc~^. 

^ Uncertainties include both statistical and system- 
atic errors : a^ = O'^^^ + Q-gy s. See Section 2.4 and 
Table 3 of |Gaztafiaga et al.| l |2009l l . 
'' Data in this row are taken from JVTS03. 
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Figure 1. Full OHD set and best-fit ACDM models. The top 
panel shows the dataset and fit results. The bottom panel shows 
the residuals with respect to the best-fit model with Gaussian prior 
on Hfj. The Hq priors are described in Scction|3.2| 



"GCH09" ) . These datasets are summarized in Table IT] 
and displayed in Figure [T] 

These datasets have seen wide application in cosmo- 
logical research. In addition to those mentioned above, 
Yi fc Zhang; (^2007j) first used the SVJ05 dataset to con- 



strain cos mological model parameters. |Samushia fc Ra^ 



tra 



( 2006 ) also used the data to constrain parameters in 

various dark energy models. Their results are in con- 
sistence with other observational data, in particular the 
SNIa. Besides parameters constraints, OHD ca n also be 
used as an auxiliary model selection criterion (Li et al. 
20091) 



Tn this paper, the OHD sets used are taken from the 
union of JVTS03, GCH09, and SJVKSIO. The SVJ05 
dataset, having been replaced by SJVKSIO, is no longer 
used, and is listed in Table [T] for reference only. Us- 
ing these datasets, we find the parameter constraints 



for a non-flat ACDM universe fi^ 



0.37 



+0.15 
-0.16 



and 



0-93-0.29 assuming the Gaussian prior Hq = 74. 2± 



3.6 km s-i Mpc-i suggested by Riess et al ( |2009| ). We 
also use a conservative, "top-hat prior on Hq, namely a 
uniform distribution in the range [50, 100], and obtained 



fin 



0.34 



+0.20 
-0.27 



and Oa = 0.86: 



-0.44 
34- 



These Hq priors 



are discussed in detail in Section [X2| and their effects on 
the parameter constraints are illustrated in Section l4] 

3. PARAMETER CONSTRAINT WITH SIMULATED 
DATASETS 

In the current absence of more OHD, we turn to simu- 
lated H{z) datasets in the attempt to explore the answer 
to the questions raised in Section [11 To proceed with our 
exploration, we must prepare ourselves with a) a way of 
generating simulated H{z) datasets, b) an "evaluation" 
model (or class of models) of cosmic expansion in which 
parameter constraint is performed using the simulated 
data, and c) a quantified measure of the datasets' ability 
of tightening the constraints in the model's parameter 
space, i.e. a well-defined "figure of merit" (FoM). These 
topics are discussed in detail in the rest of this section. 

3.1. Generation of Simulated Datasets 

Our simulated datasets are based on a spatially flat 
ACDM model with fi,„ = 0.27 and ^a = 0.73. This 
fiducial model is consistent with the 7- year Wilkinson_ 
Micro wave Aniso tropy Probe (WMA P) (Komatsu et al. 
20101), the B AO ( [Percival et al.|2010| , and SNIa (" Hicken 
et al.||20b9b ) observations. Therefore, it summarizes our 
current knowledge about the recent history of cosmic ex- 
pansion fairly well. In this fiducial model the Hubble 
parameter is expressible as a function of redshift z by 
the simple formula 



Htid{z) = HQ^n,^{i + z)^ + rtA 



(1) 



where Hq is the Hubble constant. 

The modelling of the observational data's deviations 
from the fiducial model, as well as the statistical and sys- 
tematic uncertainties of the data, can be rather difficult. 
For SNIa, there are planned projects such as the Wide- 
Field Infrared Survey Telescope (WFIRSTpl) with well- 
define d redshift distribution o f targets (|Al3cring et al 



2004 



im^ 






and uncertainty model ( |Kim et al?||2004[ [Huterer 
based on which simulated data can be generated. 
However, this is not true for OHD, as there have not 
been a formal specification of future observational goals 
known to the authors. Consequently, we have to work 



^ http : //wf irst . gsf c . nasa. gov/l 
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Figure 2. Uncertainties of H(z) in the SJVKSIO dataset. Solid 
dots and circles represents non-outliers and outliers respectively. 
Our heuristic bounds (t+ and a— are plotted as the two dotted 
lines. The dash-dotted line shows our estimated mean uncertainty 
o-Q. 



around this difHculty by approaching the problem from 
a phcnomenological point of view. 

By inspecting the uncertainties on H[z) in the 
SJVKSIO dataset (Figure I2]), we can see that the gen- 
eral trend of the errors' increasing with z despite the 
two outliers at z = 0.48 and 0.88. Excluding the out- 
liers from the dataset, we find that the uncertainties ^{z) 
are bounded by two straight lines dj^ — 16.87z -I- 10.48 
and a- = 4.41z -I- 7.25 from above and below respec- 
tively. If we believe that future observations would also 
yield data with uncertainties within the strip bounded 
by (7+ and (T_, we can take the midline of the strip 
(To = 10.64z -|- 8.86 as an estimate of the mean uncer- 
tainty of future observations. In our code, this is done 
by drawing a random number d{z) from the Gaussian 
distribution N{ao{z),e{z)) where £{z) = (ct+ — o'_)/4. 
The parameter e is chosen so that the probability of a{z) 
falling within the strip is 95.4%. 

Having found a method of generating the random un- 
certainty d{z) for a simulated datapoint, we are able to 
simulate the deviation from iJfid. Namely, we assume 
that the deviation of the simulated observational value 
from the fiducial, A7J = Hsimiz) — Ha^iz), satisfies the 
Gaussian distribution A^(0, o'(z)) from which AH can be 
drawn as a random variable. 

Thus a complete procedure of generating a simulated 
H{z) value at any given z is formed. First, the fiducial 
value iJfid(z) is calculated from equation (fll). After that, 
a random uncertainty ^{z) is drawn using the aforemen- 
tioned method. This uncertainty is in turn used to draw 
a random deviation AH from the Gaussian distribution 
A^(0, i5-(z)). The final result of this process is a datapoint 
Hsimiz) = Hfid{z) + AH with uncertainty a(z). 

In addition to the procedures described above, the 
Hubble constant (in the units of km s~^ Mpc~^) is also 
taken as a random variable and is drawn from the 



Figure 3. Snapshot of a simulated dataset realized using our 
method. The uncertainties in the simulated data are modelled 
phenomenologically after SJVKSIO, which is also shown for com- 
parison. 

Gaussian distribution A^(70.4, 1.4) suggested by 7-year 
WMAP results when we evaluate the right-hand side of 
equation (IT]). We could have fixed Hq at a constant value, 
but as we shall see in Section [3.2[ in our analysis we treat 
Hq and fJm quite differently. Namely, Q,^ is a parame- 
ter with a posterior distribution to be inferred, but Hq 
is a nuisance parameter that is marginalized over using 
some independent measurement results as prior knowl- 
edg e. This treatment c an be found in many works, such 
as ( [Stern et aL||2010a| ) and (Wei 2010]. It can be justi- 
fied oytEieneeatoreHuce the dimension of the parameter 
space given limited data, and we usually prioritize other 
parameters such as Q^ over Hq. Therefore, when gener- 
ating simulated datasets we sample Hq from a random 
distribution to reflect the uncertainty to be marginalized 
over. Otherwise, we could have "cheated" by forcing a S- 
function prior on Hq centered at its fiducial value in the 
analysis of simulated data and obtain overly optimistic 
predictions from the simulated data. 

The quality of simulated data thus generated is simi- 
lar to that of SJVKSIO. A snapshot realization of this 
simulation scheme is displayed in Figure [31 where a total 
of 128 datapoints with z evenly spaced within the range 
0.1 < z< 2.0. 

3.2. The Evaluation Model 

We use a standard non-flat ACDM model with a cur- 
vature term Jli; = 1 — 57^ — fl\ to evaluate the qualities 
of the simulated datasets. In this model, the Hubble 
parameter is given by 



H{z)^Hoy/n^{l + z)3 + fik(l + z)2 + Oa 

^HQE{z;n^,n^). (2) 

Our model choice is mainly motivated by our desire to 
reduce unnecessary distractions arising from the intrinsic 
complexity of certain cosmological models involving dark 
energy or modified gravitation. 
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We perform a standard maximal likelihood analysis us- 
ing this evaluation model. In our analysis we intend to 
marginalize the likelihood function over the Hubble pa- 
rameter Hq, thus obtaining parameter constraints in the 
(f^ni, fljy) subspace. This marginalization process also al- 
lows us to incorporate a priori knowledge about Hq into 
our analysis. 

There is a fair amount of available inform ation from 
which reasonab le priors can be constructed. ISamushia 



& Ratra (2006) used two different Gaussian priors on 
Hq: one with Hq = 73 ± 3 km s~^ Mpc~^ from 3- 



year WMAP data (S pergel et al.||2007|), the other with 
iJo = 68±4 km s-i Mpc"^ by |Gott III et al.| ( |200l| (for 
a discussion of the non-Gaus siainity of the erro r distribu- 
tion in Hq measurements, see Chen et al. 2003 ). Lin et al. 



( 2009J usedgp == 72± 8 km s"^ Mpc"^ as suggested by 
( I'reedman et al.|2001 ) . In our work we use a more recent 
determination of Hq = 74.2 ± 3.6 km s~^ Mpc~^ (Riess 
et al. 2009 1 as an update to the ones cited above. We also 
consider a "top-hat" prior, i.e. a uniform distribution in 
the interval [50, 100]. Compared with any of the peaked 
priors above, this prior shows less preference of a particu- 
lar central value while still characterizing our belief that 
any value of Hq outside the said range is unlikely to be 
true. Notice that the intrinsic spread of Hq involved in 
the generation of simulated datasets is smaller than ei- 
ther prior chosen in this step, in consistent with the belief 
that the prior adopted in the estimation of parameters 
should not be spuriously optimistic. 

Having chosen the priors on Hq, it is straightforward to 
derive, up to a non-essential multiplicative constant, the 
posterior probability density function (PDF) of parame- 
ters given the dataset {Hi} by means of Bayes' theorem: 



P{n^, n^l {H,}) = J P(a„, ^A, Ho\ {H,}) diJo 

= J C ({iJjIf^n., f^A, Hq) P{Hq) dHQ 



(3) 



where C is the likelihood and /'(-ffo) is the prior PDF 
of Hq. Assuming that each measurement in {H^} has 
Gaussian error distribution of H{z) and is independent 
from other measurements, the likelihood is then given by 



c{{H,}\n^,n^,Ho)^ (H 



1 



where the x^ statistic is defined by 

2 _ v^ [HQE{zi; rim, Q\) - 



cxp 



H,: 



XT 
' 2 



(4) 



(5) 



and CTi's are the uncertainties quoted from the dataset. 
The posterior PDF can thus be found by inserting equa- 
tion Q and the exact form of P{Hq) into equation I^ 
We now show that the integral over Hq in equation (p| 
can be worked out analytically for our P{Hq) choices. 



Gaussian prior. 

P{Hq) : 



Let 



v/2^ 



exp 



{Hq - hh) 
2a\ 



(6) 



In this case, equation (pi reduces to 

p(r!„„f7A|{-ffj) 



erf 



B 



exp 



52 

A 



(7) 



where 



1 



._ ^ M-^ E'^{zi;n,-n,nA) 

E{zf,n^,nA)H, 



B = 



1a\ 



+E 



la 



and erf stands for the error function. The form shown 
in equation ^\ is not normalized; all multiplicative con- 
stants have been discarded. 

Top-hat prior. — The PDF of a uniform distribution over 
the interval [x, y] can be written as P{Hq) = Q{Hq — 
x)Q{y — HQ)/{y — x) where 6 denotes the Heaviside unit 
step function. With this prior on Hq, equation ([3| be- 
comes 



_U{x,C,D)-U{y,C,D) 



where 
and 



E 



E''(z,-n^,nA) 

2a? 



G 



D 



exp 



~C 



(8) 



E 



E{zi;Q.m,Q.A)Hi 



2a 



U{x,a,f3) — erf 



/3 — xa 



la 



The normalization constant has been dropped from for- 
mula ([8| as well. 

3.3. Figure of Merit 

The posterior probability density functions obtained 
above put statistical constraints over the parameters. As 
the dataset {Hi} improves in size and quality, the con- 
straints are tightened. To evaluate a dataset's ability 
of tightening the constrains, a quantified figure of merit 
(FoM) must be established. 

We note that the FoM can be defined arbitrarily as 
long as it reasonably rewards a tight fit while punishing 
a loose one. Its definition can be motivated purely sta- 
tistically, for example the reciprocal hypervolu me of the 
95% confide nce region in the parameter space (Albrecht 
et al. 20061. However, a definition that is sensitive to 
some physically significant structuring of the parameter 
space can be preferable if our scientific goal requires it 



of Albrecht et al. 



( |Linder)[2006l ). 

In this paper we use a st atistical FoM de fi nition sim- 
ilar to the ones 
(2008T), and ^ 
is defined to 

p(r!„,i7A|{a,}) 



|Bueno 



S anchez et ai. 
enciosec 



the area 

exp(-AxV2)P, 



2006 ), 
(12^. 



Liu 



et al. 
F5M 



The contour of 
where Pmax is 
the maximum value of the posterior PDF, and the con- 
stant Ax^ is taken to be 6.17. This value is so chosen that 
the region enclosed by this contour coincides with the 2cr 
or 95.4% confidence region if the posterior is Gaussian. 
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In the rest of this paper we will use the term "2(T region" 
to refer to the region defined in this way. The la and 3a 
regions can be defined in the same manner by setting the 
constant Ax^ to 2.3 and 11.8 respectively. It is important 
to bear in mind that our definition of "na regions" are 
motivated by nomenclature brevity rather than mathe- 
matical concreteness, since the posterior PDFs (eqs. [t] 
and Is]) are manifestly non-Gaussian. 

We make two further remarks on our FoM definition. 
First, since our definition of FoM is statistical rather than 
physical, we do not exclude the unrealistic parts of the 
confidence regions from the area. Second, the FoM by our 
definition is obviously invariant under the multiplication 
of P(nm, ^a| {Hi}) by a positive constant, which ju stifie s 
the omission of normalization constants in Section 13.21 



4. RESULTS FROM THE SIMULATED DATASETS 

Using the method described in Section [XT] we gener- 
ate 500 realizations of the simulated H{z) daiaset. Each 
realization contains 128 datapoints evenly spaced in the 
redshift range 0.1 < z < 2.0. From each realization, 
successively shrinking subsamples of 64, 32, and 16 dat- 
apoints are randomly drawn. These subsamples are then 
used in conjunction with the full OHD set to obtain their 
respective figures of merit. 

The figures of merit are naturally divided into 4 groups 
by the size of the simulated subsample used in the cal- 
culation. For each group, median and median absolute 
deviation (MAD) statistics are calculated. The median 
and the MAD are used to represent the central value 
and spread of FoM data respectively, and they are used 
in preference to the customary pair of the mean and 
the standard deviation, because they are less affected 



by egregious outliers (see NIST 2003 Chapter 1.3.5.6). 
The outliers mainly arise irom worst case" realizations 
of simulated H(z) data, for example one with a large 
amount of datapoints deviating too much (or too little) 
from the fiducial model. 

To compare the parameter constraining abilities of our 
simulated H{z) datasets with those of SNIa data, we 
have also fitted our evaluation model (eq. [2]) to the 
ConstitutionT dataset which is a subset of theCons titu- 
tion redshift-distance dataset ( [Hicken et al. 2009b) de- 
prive d of outliers that account for internal tensions ( Wei 
20101. It is worthwhile to point out that the prior or Hq 
used in the SNIa fitting procedure is fundamentally dif- 
ferent from either one discussed in Section [X2| Namely, 
when SNIa data is used, the parameter Hq and the in- 
trinsic absolute magnitude of SNIa, A/q, are combined 
into one "nuisance parameter" M = Mq — 51g-ffo, and 
is marginalized o ver under the as sumption of a fiat prior 
over (— oo,-|-cx)) (jPerlmutter et al. 1999). This discrep- 
ancy should be kept in mind when comparing or com- 
bining SNIa and H{z) datasets, and we hope it could be 
closed in the future by better constraints, cither theoret- 
ical or observational, on Mq. 

Our main results are show in Figure ]4] As one may 
intuitively assume, the median FoM increases with the 
size of the dataset. We find that the data subset with 64 
simulated datapoints leads to a FoM of 8.6 ± 0.7 under 
the Gaussian prior on Hq. This median FoM already 
surpasses that of ConstitutionT. However, the top-hat 
prior leads to significantly lower FoM. Under the top-hat 
prior we used, as many as 128 simulated datapoints are 
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Figure 4. Figures of merit from each group found by combining 
all 500 realizations. FoM medians are plotted against the sizes 
of the simulated subsample, and median absolute deviations are 
shown as error bars. Horizontal lines across the figure marks the 
figures of merit of purely observational datasets. 
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Figure 5. Confidence regions in the (f2ni,f2A.) parameter sub- 
space calculated using the snapshot realization shown in Figure Is] 
The Gaussian prior Hq = 74.2 it 3.6 km s~^ Mpc~^ is assumed! 
The shaded regions, from lighter to darker, correspond to the 2cr 
constraints obtained from the data subsets containing 16, 32, 64, 
and 128 simulated H{z) datapoints respectively. For comparison, 
the 2(7 region from pure OHD alone, with the same Gaussian prior 
on Hq, is plotted in the dashed contour. The three dotted con- 
tours are the 1, 2, and 3cr constraints from the ConstitutionT SNIa 
dataset. The solid straight line signifies the boundary of fii^ = 0. 



needed to reaches the FoM level of ConstitutionT. 

The degeneracy of confidence regions obtained from 
H(z) data in the (fim, ^a) parameter subspace is shown 
in Figure [5] Because of this degeneracy, H{z) datasets 
cannot be used alone to constrain ilk well. This degen- 
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erate behavior is similar to that of SNIa dat£ 



5. FISHER MATRIX ANALYSIS OF FUTURE DATA 

The simulated data used in Section [3] are based on 
the quality of currently available measurements, and we 
have tried not to be too optimistic about their uncertain- 
ties. However, we have reasons to expe ct an increase in 
the qu ality of future H{z) data. First, Crawford et al 
(2010) analysed the observational requirement of mea- 
surmg H{z) to 3% at intermediate redshifts with age- 
dating. Second, the Baryon Oscillation Spectroscopic 
Survey (BOS^j) is designed to constrain H{z) with 2% 
precision at redshifts z w 0.3 and 0.6 by measuring BAO 
imprints in the galaxy field, and at z « 2.5 using the 
Lyman-a absorption spectra of quasars. 

By incorporating these specifications of future data, 
we can estimate their expected fig ures of merit using th e 
Fisher matrix forecast technique (Dodelson et al. 1997 ). 
The 3x3 Fisher matrix F is calculated from equation ^ 
with the (Ti's determined by future data specifications, 
and the evaluation of matrix elemen ts is made at the 
fiducial parameter values (see Section 3.1). The matrix 
elements are the second partial derivatives of x^ with 
respect of the parameters: 



1 d^x' 



^'^ 2 ae.de, 



(9) 



where the S^'s are the parameters, namely (r^m, Q\, Hq). 
Notice that we use the word "Fisher matrix" only loosely. 
The second derivative matrix defined by equation Q is 
not the Fisher matrix in t he strict sense, but the laeas 
are intimately related (see Dodelson||2003 ). 

In order to obtain the l-olVl in the 2-dimensional pa- 
rameter space of (riin,^A)j we must marginalize over 
Hq. We adopt the Gaussian prior on Hq with an = 
3.6 km s~^ Mpc~^, which is the same as the one used in 
Section |3.2 Gaussian marginalization is performed us- 
ing a straightf orward m odification o f the projection t ech- 
nique from ( [Dodels on "2003) and ( [Press e t al. '2007'). A 
more ge neral and detailed analysi s of this problem was 
made by Taylor & Kitching (2010), but for our purpose, 
simply adding l/c^ to F33 and then projecting onto the 
first two dimensions will do the work. This is applicable 
when the mean of the prior Hq is close to the fiducial 
value, and we have numerically verified the validity of 
this approximation in our work. 

Denoting the marginalized Fisher matrix by F, the iso- 
Ax^ contour in the parameter space is approximated by 
the quadratic equation 



{AOfFAe = Ax^ 



(10) 



where A6 = 6 6f^^ is the parameters' deviation from 
the fiduci al va lue, and Ax^ is our chosen constant 6.17 
in Section 3.3 By construction, F is positive-definitq^ 
therefore the above equation describes an ellipse. Its 



enclosed area is simply A — tt / y dei{F / Ax^) , therefore 



* http : //www . sdss3 . org/boss . php 

^ I'his can be proved u sing the pos itivc-dcfiniteness conditions 
of the Schur complement (IPuntanen &; Styan 2005). 
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Figu re 6. Predicted FoM using the error model of jCrawford et al.| 
|2010[ |. The dependence of FoM on the number of nieasurements is 
shown for three cases: the solid, dashed, and dotted lines are for 
3%, 5%, and 10% relative error respectively. The horizontal line 
across the figure shows the FoM of ConsitutionT for comparison. 



we can estimate the FoM by 



^°^ = i = y^^KA^- 



In our setup, we assume that t he rel ative error of H{z) 
anticipated by Crawford et al. (2010) could be globally 
achieved within the redshift interval 0.1 < z < 1.0. Un- 
der this assumption, wc can estimate how many data- 
points will be needed to reach the ConstitutionT FoM 
using the Fisher matrix method discussed above. We 
chose not to incorporate the available data so that we 
can work with the simple error model specified by future 
data. 

Our main results are shown in Figures [6] and [T] For 
relative H{z) errors of 3%, 5%, and 10%, the required 
number of measurements are 21, 62, and 256 respectively. 
This number N increases steeply as the relative error 
increases, as shown in Figure [7] 

We can also apply the Fisher matrix analysis to future 
BOSS data. We find BOSS alone gives FoM « 15. This 
promising result is a combined consequence of its high 
precision and extended redshift coverage. In Figure [8] we 
plot the Fisher matrix forecast of the confidence regions. 



6. CONCLUSION AND DISCUSSIONS 

We have explored the possibility of using OHD as an 
alternative to SNIa redshift-distance data in the sense 
of offering comparable or higher FoM. By using simu- 
lated H{z) datasets with an empirical error model sim- 
ilar to that of current age-dating data, we show that 
more than 60 future measurements of H(z) in the redshift 
range < z < 2 could be needed to acquire parameter 
constraints comparable with those obtained from SNIa 
datasets like ConstitutionT. In addition, precise and ac- 
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Our result furthers a conclusion of Lin et al. (20091 

namely, that OHD play al- 
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Figure 7. Number of measurements required to matcfi Constitu- 
tionT, TV, as a function of the relative error of H{z). The curve 
turns upward steeply when the error is large. 
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Figure 8. Same as Figure!?] but for the forecast of future data. 
The light shaded region is the la region from the 21 ago-dating 
measurements with relative error 3%. The unfilled region bounded 
by the dashed contour is the 2a region from BOSS H{z). The 
joint 2a region of BOSS and age-dating is shown as the inner and 
darker shaded area. Like in Figure!?] 1, 2, and 3a constraints from 
ConstitutionT are shown as dotteocontours for comparison. 



curate determination of Hq is crucial for improving the 
FoM obtained from OHD, and a broad prior on Hq leads 
to lower FoM. When we progressively lower th e error 
of future measu rements to 3% as discussed by (Craw- 
ford et al.|2010 l, we estimate that ~60 measurements m 
the shorter redshift interval [0.1,1.0] will be needed to 
achieve the same result. In summary, we give an affir- 
mative answer to the first of the two questions raised in 
Section [l] and a semi-quantitative answer to the second. 



and Carvalho et al. 



02008] ) 

most the same role as that of SNIa for the joint con 
straints on the ACDM model. We have shown that the 
OHD set alone is potentially capable to be used in place 
of current SNIa datasets if it is large enough. We note 
that our forecast of OHD data requirement is competitive 
in the sense of observational cost compared with super- 
nova observations. Throughout our analysis, we used the 
ConstitutionT dataset used as a standard of FoM. This 
dataset is a subset of the Constitution compilation, a 
combina tion of the ground-b ased CfA3 supernova obser- 
vations ( Hicken et al.]]2009a ) and Union, a larger compi- 
lation of legacy supern ovae and space-based observations 
( Kowalski et al.][2008 |. The CfA3 sample alone requires 
: nights tor each ot the 185 supernovae observed. On 
the other hand, current OHD from age-dating does not 
require space-based observ ations. As described in de- 
tail by Stern et al. (2010b I, 24 galaxy cluster containing 
target chronometer galaxies were obtai ned in only two 
night s using the Keck I telescope. In (Crawford et al. 
20101, it was estimated that the South African Large 
I'elescope (SALT) is capable to measure H{z) to 3% at 
an individual redshift in ~180 hours. 

Admittedly, our approach to the problem is tenta- 
tive as well as model-dependent. The conclusion is 
reached under the assumption of the fiducial ACDM 
model (e q. fll ), the model of uncertainties described in 
Section s |3.1] ~and J5] as well as the evaluation model of 
Section 13.2] Therefore, care must be taken not to ex- 
trapolate our conclusions well beyond those assumptions, 
especially when used in the planning of observations. 

Even if our uncertainty models fit real-world observa- 
tions well, in practice the results may be less promising 
than what we suggested in this paper. We assumed a 
fairly deep redshift coverage which may be difficult to 
reach by current observations (see ]Stern et al.]]2010b 



for the redshift distribution of LRCs used to deduce 
SJVKSIO). Ultimately, the FoM from OHD may be lim- 
ited by our ability to select enough samples for future 
observation and the redshift distribution of these sam- 
ples, rather than by the relative error of each individual 
measurement. Fortunately, the BOSS project could ex- 
tend H{z) measurement into deeper redshift. Also, It is 
worth notin g that the proposed Sandage-Loeb observa- 



tional plan ( Corasaniti ^et al.j 2007 Zhang et al. 2010 
Araujo fc Stoeger^ 2010' ) could be used to exteiid our 



knowledge ot cosmic expansion into the even deeper red- 
shift realm of 2 < z < 5 by measuring the secular varia- 
tion of cosmological redshifts of the Lyman-a forest with 
future high-resolution, space-based spectroscopic instru- 
ments (see (Sandage 1962) and (Loeb 1998) for the foun- 
dations of the test). 

Finally, we note that future CMB observation pro- 
grams, such as the Atacama Cosmology Telescope]^ may 
be able to identify more than 2000 passively evolving 
galaxies up to z ~ 1.5 via the Sunyaev-Zel'dovich ef- 
fect, and their spectra can be analyzed to yield age mea- 
surements that will yield app roximately 1000 H (z) de- 
terminations with 15% error (Simon et al. 2005| ). This 
promises a future data capacity an order ot magnitude 
more than what we have estimated to be enough to match 



http://wMW.physics.princeton.edu/act/index.html] 
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current SNIa datasets. Combining this prospect with 
future high-z, high-accuracy H(z) determinations from 
BAO observations, it is reasonable to expect that the 
OHD win play an increasingly important role in the fu- 
ture study of the expansion history of the universe and 
cosmological parameters. 
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